Impurity effects on semiconductor quantum bits in coupled quantum dots 
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We theoretically consider the effects of having unintentional charged impurities in laterally coupled 
two-dimensional double (GaAs) quantum dot systems, where each dot contains one or two electrons 
and a single charged impurity. Using molecular orbital and configuration interaction method, we 
calculate the effect of the impurity on the 2-electron energy spectrum of each individual dot as well 
as on the spectrum of the coupled-double-dot 2-electron system. We find that the singlet-triplet 
exchange splitting between the two lowest energy states, both for the individual dots and the coupled 
dot system, depends sensitively on the location of the impurity and its coupling strength (i.e. the 
effective charge). A strong electron- impurity coupling breaks down equality of the two doubly- 
occupied singlets in the left and the right dot leading to a mixing between different spin singlets. 
As a result, the maximally entangled qubit states are no longer fully obtained in zero magnetic field 
case. Moreover, a repulsive impurity results in a triplet-singlet transition as the impurity effective 
charge increases or/and the impurity position changes. We comment on the impurity effect in spin 
qubit operations in the double dot system based on our numerical results. 

PACS numbers: 73.21. La, 03.67.Lx, 73.23.Hk 



I. INTRODUCTION 



The goal of this work is the calculation of the low- 
lying energy spectra of 2-electron semiconductor quan- 
tum dots (QDs) (both single dots and coupled double- 
dots) in the presence of nearby static charged impurity 
centers within a minimal model. The purpose is to quan- 
tify the effects of quenched random charged impurities 
on the singlet-triplet splitting in QDs in order to assess 
the importance of background unintentional impurities, 
which are invariably present in the environment, in ad- 
versely affecting the operations of exchange-coupled dots 
as elementary spin qubits for solid state quantum com- 
putation. Since the positions and the strengths of the 
background unintentional impurities are unknown and 
random, we study the impurity effects as functions of 
impurity position and coupling strength (defined as the 
effective impurity charge Z) assuming the impurities to 
be Coulombic charge centers so that their effective po- 
tential falls off slowly as a 1/r potential away from the 
impurity location where 'r' is the distance from the im- 
purity. Since the main background impurities in GaAs 
and Si, the two most relevant semiconductors of interest 
for solid state spin quantum computation, are random 
charge centers, our consideration of Coulombic impurities 
with a long-range impurity potential is reasonable. The 
theoretical results presented herein, while being com- 
pletely microscopic and fully quantum mechanical, are 
phenomenological in nature since the impurity position 
and charge are treated as unknown parameters. While 
our results show clearly the strong effect of local back- 
ground charged impurities on the low-lying QD spectra, 
its usefulness is limited in comparing with experiment 
since no direct information about impurity locations is 
currently available experimentally. On the other hand, 
our results establish the manifest importance of using the 
highest quality background materials for semiconductor 



spin qubit operations since the presence of even a sin- 
gle charged impurity in the vicinity of the QDs seems 
to completely ruin the operational logistics of coupled 
QD systems. The exchange coupling (or equivalently, the 
singlet-triplet splitting) depends strongly and sensitively 
on the location and the strength of the charged impurity, 
which means that (1) a single charged impurity located 
nearby could destroy the qubit, and (2) even a remote 
charged impurity could have a strong adverse effect in- 
ducing substantial qubit decoherence if the impurity is 
dynamic and has a fluctuation timescale comparable to 
gate operations timescales- in fact, the impurity fluctua- 
tion timescale will become a dominant decoherence time 
since the exchange energy will vary substantially over 
this timescale. The motivation of our work is a clear un- 
derstanding of the energetics of QD systems in the few 
electron situation in the presence of charged impurities 
so that some rudimentary quantitative magnitudes of the 
impurity effects are available for qubit operation consid- 
erations. 

Coupled QDs for quantum computationi^-S have been 
extensively studied due to the prospect of using QDs as 
ideal environment to confine and manipulate the QD 
electron spins. The quantum bit, or qubit, of infor- 
mation is encoded and stored in these localized single 
electron spins which exploit a spin relaxation time of 
the order of milliseconds^^-, sufficiently long to allow 
the performance of coherent spin operations. The pro- 
posed quantum computer^ in solid states operates based 
on the exchange coupling J between the two electron 
spin qubits manipulated by an external magnetic field. 
This exchange energy can be envisioned as the effective 
coupling between the two spins in the double dots^ via 
the Heisenberg spin Hamiltonian, H ~ Jsi • S2, which 
takes into account possible contributions from different 
hybridzation between singlet and triplet states. Hence, 
J is determined through the gate voltage control over 
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the tunneling coupling between the coupled two QDs. A 
complete understanding of J is important because it di- 
rectly determines the a/ SWAP operator which describes 
the exchange information between the two qubits in the 
double dots. The fast solid-state two-qubit operation''^, 
generated as a consequence of the electron spin exchange 
under electrical control, and the single-qubit operationl*^ 
suffice to assemble a standard quantum computing sys- 
tem. The number of electrons in such a QD system can 
be controlled preciselj*ii"— to 0, and the electron ex- 
change interaction is tunable by the applied gate volt- 
ages. Thereby, the coupling between the dots can also be 
controlled. 

Established quantum regimes such as quantum entan- 
glement between individual electrons in one dot with the 
other electrons in the other dots and superposition of 
electron spin qubits are the major objects^i^i^ii^"— in 
such exchange coupled QD system. Prerequisite crite- 
ria for realization of a quantum computation system'^ 
(including initialization, manipulation of spin qubits, 
and readout) have been demonstrated for single dot^^ 
and coupled QD system^-, provided the long relaxation 
time^ii^ of the spins, by using charge sensing and fast 
spin-to-charge conversion techniques. One of the per- 
ceived advantages of solid-state quantum computing is 
the scalability with existing semiconductors. In ad- 
dition, integrating multi coupled dots into a quantum 
circuit is made possible by adding more suitable gate 
electrodes fSilS, 

QDs in versatile GaAs semiconductors not only have 
been considered the most widely-studied objects in the 
QD science but also their well-understood physics are ap- 
plicable to a variety of materialsj^i Often, unintentional 
impurities found in the dot sample are used in the fabrica- 
tion process to adjust the potential well height between 
different heterostructure semiconductors (such as Si in 
gated GaAs/AlGaAs QDs - see e.g. Refi^) allowing the 
charge flow of electrons. Statistically, such QD system 
containing impurities can be studied for one, two, ... im- 
purities. On the other hand, when integrating multiple 
coupled 2-dot systems utilized as a multi-qubit gate, in 
the non-overlapping regime between different 2-dot sys- 
tems, a single spin qubit in one coupled 2-dot system can 
be affected^^ by the other coupled 2-dot systems. Each 
coupled 2-dot system acts as a source of electrostatic field 
to the others and can be treated as "charged impurities" . 

Theoretical studies of impurities in coupled QDs have 
been rarely found in the literature. In fact, the impurities 
are practically found randomly in/outside the dot sam- 
ple and their positions cannot be specified precisely. In a 
coupled triple-dot system, a relatively large collection of 
impurities was considered^^. These statistical impurities 
were theoretically assumed to induce a weak perturba- 
tion to the coupled triple QDs. The authors^^ found that 
the Coulomb exchange energy of the impurities with the 
QD electrons in this study resulted in decoherence of the 
coded qubit states. Thus, any information processing us- 
ing electron charge degree of freedom needs to take into 



account the decoherence channel due to charge fluctua- 
tions. 

In the present paper, we study the influence of charged 
impurities on the singlet-triplet splitting of the two lowest 
energy levels as well as the energy spectrum as a function 
of the impurity position, the impurity effective charge, 
and the number of impurities for a coupled two-dot sys- 
tem in zero magnetic field. We also examine the impurity 
effect on the coupled electron qubits by tuning the po- 
tential well height to different values and obtain different 
triplet-singlet transitions for the repulsive impurity case. 
To accomplish, we discuss the singlet-triplet splitting in 
the presence of two impurities with similar charge located 
in the two separate dots of the system. 

The paper is organized as follows. In Sec. II we in- 
troduce the model and methodology. Impurity effect on 
energy spectrum of a two-electron single QD containing a 
single charged impurity is first discussed in detail in Sec. 
III. Sect. IV is spent to present our studies on a cou- 
pled 2-dot system in the presence of one or two charged 
impurities. We examine the impurity position and impu- 
rity effective charge dependence of the singlet-triplet spin 
splitting energy. Influence of the confining potential bar- 
rier height on the electron-impurity (e-I) and electron- 
electron (e-e) coupling is also explored. All the results 
presented in this paper are obtained at zero magnetic 
field, B—0. Summary of our results and conclusion are 
found in Sec. V. 



II. THEORETICAL MODEL 

Horizontally-coupled QDs are grown by the depletion 
of the two-dimensional electron gas (2DEG) using typ- 
ically the gated mechanism^i^. Such gated QDs have 
typical size of about few tens of nm. Consequently, the 
lowest excitation energy in such a QD system is found 
of the order of few meV. When the inter-dot coupling 
strength is substantial, the electrons in the coupled dots 
strongly quantum mechanically couple with each other. 
As long as the phase is coherent, the electrons can "tun- 
nel" between the adjacent dots forming different entan- 
gled qubits. Using single-conduction-band effective-mass 
approximation which was justified in Ref.^^, the Haniilto- 
nian describing a two-electron coupled double QD system 
containing unintentional charged impurities (charge size 
Ze) in zero magnetic field reads: 



H = hi+h2 + Vc + K 



e-Imp 



I-T 



(1) 



Here, the first two terms hij are the single-particle 
Hamiltonian of the two quantum-dot electrons (coordi- 
nates Ti) 



2m* 



(2) 



confined by a potential well V{r). In the present pa- 
per, this confining potential is constructed as a linear 
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combination of the three different Gaussians and can be 
separated into two parts: 



V{v) = Va{e 



~P ^1^' , ' — '^7 



V) (3) 



The first part acts as a double-well confining potential 
for the coupled double dot system and the second part is 
used to control the electrostatic barrier height between 
the two dots independently. The set of varying param- 
eters Vo,Vb,lxily,hx,hy characterizes the potential well 
depth and barrier height. A nonvanishing overlap be- 
tween the wave functions of the two quantum dots signi- 
fies the electron virtual tunneling between the two dots, 
i.e. the exchange energy is nonzero. Because this prop- 
erty can be tuned by the applied gate voltage this ex- 
change energy thereby can be controlled. It is worth 
noting that Vt can independently modulate the barrier 
height resulting in concomitant change in the inter-dot 
separation without modifying the single-particle energy 
spectrum of the individual QDs in the system. The last 
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FIG. 1: (Color online) Confinement potential double well of 
coupled double QDs modeled along the a;-direction plotted for 
dilTerent barrier depth 14 (from 35 down to 11.4 meV) (left 
column). The corresponding double well minima (±0^,0) 
are indicated in the same line in the right column of the ref- 
erence. The confinement energy at the two well bottoms Em 
is enclosed as the magenta open circles for various values of 

three terms in Eq. ([T]) are the Coulomb interaction be- 
tween, respectively, the two electrons 

Vc 



47reeo|rij| 
the electron and the impurities 

lS!i=2 N^=2 



e-Imp ^ 4.e.o|R.~n| 



(4) 



(5) 



where Ni and denote the number of impurities and 
electrons, respectively, and the impurities 



Z1Z2 



l-l 



47reeo|Ri -R2I 



(6) 



The solution of a single particle confined by a parabolic 
potential is well known as the Fock-Darwin basis: 



1 



with corresponding energy 

E^,i = hLOo[2n+\l\ -t- 1) 



(7) 



(8) 



where [n, I) stand for radial and azimuthal quantum 
number, respectively. Iq is the confinement length which 
is defined via the confinement frequency ujq: Iq = 

— r — . Using this basis as the radial part we can con- 

struct the many-body wave function for the considered 
system. 

We introduce a dimensionless parameter, 



(9) 



with at, ~ the effective Bohr radius, which is used 

to describe the relation between the effective Rydberg 

* 4 

energy R* — 2h'Hi7veoey^ ^^'^ ^^'^ confinement energy Hluq 



(10) 



Both e-e and e-I Coulomb interaction are evaluated in 
terms of matrix elements as 



(*l^e-Impl*> = ^o''(*|- 



1 



I*) 



I*) 



(11) 



where ^I* denotes the wave function of the system, r^j — 

r/lo, R = R//o, and Vq'' — e^/47reeo^o is the Coulomb 
energy unit. V(f relates to the confinement energy via 
the relation: 



A • huji 



0- 



(12) 



The numerical results are implemented for e.g. a GaAs 
QD with TO* = 0.067mo, e = 13.1, ge = -0.44, i?,* = 5.31 
meV (corresponding to a*Q ~ 10.3 nm). A is changed 
by changing the confinement energy hioo. A = 1 gives 
huJo 

length ^0 — o-B 

For Si/SiGe (e=13) and Si/SiOs (e=6.8) quantum dots 
with a heavier effective mass m* = 0.19too, effective 



2i?* = 10.62 meV and corresponding confinement 
10.3 nm. 
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Bohr radius and effective Rydberg energy are, respec- 
tively, = 4 nm and R* — 15 meV, and a*g = 2.11 
nm and R* = 44.76 meV (see e.g. Ref.— ). A = 1 gives 
hujo = '2R; = 30 meV {Iq = 4 nm) for a Si/SiGe QD and 
hojo = 2i?; « 90 meV (Iq = 2.11 nm) for a Si/SiOs QD 
system. Even though our numerical results are applica- 
ble to GaAs QD systems, singlet-triplet splitting energies 
for Si QD systems are also provided in Appendix |D] for 
comparison purposes. 

Reducing A means that the effective length, Iq, of the 
QD decreases while the energy spacing, huo, between the 
2D shells, i.e. the s-, p-, d- levels, will increase. In the 
small A limit, the problem at hand converts to the prob- 
lem of independent particles. In the opposite case, very 
large A, the problem approaches the classical situation. 

In the coupled double QD system, the single-particle 
solutions in each dot are obtained approximately based 
on an assumption that around the center of the each 
dot (±am,0) the single-electron problem can be treated 
as a 2D harmonic oscillator. This means that the con- 
fining potential well V{r) performs a quadratic form 
V{r) — Em « ]^[(a; ± a™)^ -I- y^], with Em the bottom 

energy of the potential well, around (±am,0). Changes 
in Em when Vb varies can be obtained as the open cir- 
cles in Fig. [T] The single-particle wave functions now are 
identical to the Fock-Darwin levels centered at (±am, 0) 
and the single-particle energy spectrum is shifted by an 
amount of Em ■ 



1 



^L{R) (r) = j- 



iW 



lo V 7r(n+ |Z|)! V / 



(13) 



of our calculations except when we examine the barrier 
height dependence of the exchange energy J where Vf, 
can vary. The center region of each QD, however, has 
unchanged effective length, namely 1^^ regardless to the 
change in the barrier height. The other parameters taken 
after Refsi^ '^^i^^ are Ix — ly — 30 nm, a ^ Ih^ — Ix/^ — 15 
nm, Iby = 80 nm. 

We assume that the impurities are located arbitrarily 
in or outside the coupled QDs. Their coordinates are 
R-fe = ixk,yk, Zk), {k=l,2}. Theoretically, the effective 
coupling between electron and the localized impurities 
as well as the coupling between the impurities with each 
other can be tuned by engineering the impurity charge 
Z. 

Configuration interaction (CI) and molecular orbital 
(MO) method are used to numerically solve the Hamilto- 
nian Eq. ([T]). Both construct the total wave function 
of the system as a superposition of different possible 
quantum configurations (CI) or molecule states (MO) ex- 
tended in the basis of single-particle wave functions: 



*(ri,r2) = ^Vz(ri,r2) 



(15) 



where V'i(ri,r2) represents one many-electron configura- 
tion as a Slater determinant. Each term of this Slater 
determinant is a single-electron wave function consisted 
of the radial part, the Fock-Darwin state </?(ri,r2) defined 
in Eq. (|13p . and the electron spin part (detailed justifi- 
cations are referred to Refi^). 

The singlet-triplet spin splitting energy of the electrons 
J in such a coupled 2-dot system is defined as the energy 
difference between the two lowest singlet {'^^) and triplet 
(^f^) states: 



where tl,r. = {x ± am,y) and ojq = y 
quadratic confining frequency which defines a new length 
?n = —T- — called the confinement length. The single- 
particle ground-state wave function is 



J = 



Hund-Mulliken 



(16) 



[(rr±a,„)^ + B^l 



(14) 



To quantitatively evaluate the advantage of using the 
above confining potential present in Eq. ([3]), we plot 
in Fig. [T] the confining potential well where its barrier 
height and the QD centers are modified by changing Vt- 
Here, Vb is reduced from 35 meV to 30, 25, 20, 15, and 
11.4 meV. As a result, the barrier height will decrease 
making the electron exchange energy increased. For ex- 
ample, the system with Vb — 35 meV has the corre- 
sponding Gm ~ 17.52 nm (« 1.75^o), and barrier height 
AVfc — 13.26 meV. For Vb = 30meV, these parameters 
are a™ « 16.41 nm (« 1.64^0) and AVb = 9.65 meV. 
Decreasing Vb leads to a shortened inter-dot separation 
2am and a smaller AVJ,. Details can be obtained in Fig.[T] 
In our numerical results, we use Vb = 30 meV for most 



In Hund-Mulliken model, the energy spectrum consists 
of four levels which arc four possible superpositions of the 
four basis (three singlets and one triplet) wave functions 

V'f(i-i,r2) = -^[(^L(ri)(y3K(i-2) + 'PL(ri)(/3fl(r2)](17) 

V'|'(ri,r2) = tpLiri)ipL{r2) 
V'f(ri,r2) = ipR{ri)tpR{r2) 

V'^(ri,r2) = ■^[(^L(ri)v3fl(r2) - (pL(ri)v?fl(r2)]. 

In the coupled QD system without an impurity, these 
three singlets do not couple with the maximally entangled 
triplet state ip'^ therefore they can be treated separately. 
The entire Hamiltonian matrix represents these singlets 
and triplets as independent blocks. 
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B. Configuration Interaction 

The most difficult task in finding the eigenvalues of 
a coupled 2-dot system lies in the basis choice among 
which Fock-Darwin basis and Gaussian basis are the most 
widely used. However, in both cases a closed analytical 
form for the Hamiltonian matrix elements, essentially the 
e-e Coulomb matrix elements, has not yet been obtained. 
The reason is that the single-particle solutions of differ- 
ent dots i) have distinct zero-points shifted to the two 
bottoms of the confining potential well and ii) thus are 
not orthogonal with each other. Consequently, the num- 
ber of distinguishable single-particle quantum states for 
the double-dot system will be doubled. Effectively, the 
size of the entire Hamiltonian increases in comparison to 
the single dot case. Quantitatively, for 3^ = subspace, 
such a number is 4 times larger than that of the single 
dot problem (2NX2N). Specifically, if only the s-waves 
are taken into account, the number of configurations in 
the subspace = is 1 for single QD and 4 for coupled 
2-dot system and if the s- and p-waves are included those 
numbers are 9 and 36, respectively. 

The above fact poses much difficulty to solving the 



eigenvalues of the coupled QD Hamiltonian. The most 
time consuming part is spent to calculate the Coulomb 
matrix elements. In the problem at hand, the number 
of Coulomb elements increases due to the e-I exchange 
interaction. However, based on the following facts one 
can take out nonphysical excited single-particle quantum 
states: i) the form of the confining potential, which is 
tuned electrostatically by metal top gates, is not exactly 
known and ii) the barrier potential height is finite. The 
latter factor validates the harmonic approximation for 
the confinement potential double well resulting in only a 
limited number of Fock-Darwin states involved. 



III. IMPURITY EFFECTS IN SINGLE 
QUANTUM DOTS 

In the single QD case, we assume that there is only 
one impurity and the impurity is located along z-axis, 
i.e. R = (0,0, d). The addition Coulomb interaction of 
the electrons with the impurity is obtained analytically 
(see Appendix [M : 



2 r(ni+/+-n)r(n2+/+ + i) 

A y^niW-ini + Z+)!(n2 + /+ + 1)! 



(18) 



V- V (-^i). {-n2)kr{l+ +j + k + l)r(l+ + 1) ~ 

jikir{i+ + j + i)rii+ + k + i) ""^ ' 



j=0 fe=0 



where Z+ = \li\ = IZ2I and 

Integral (|19p can be obtained through the recurrence: 
Io{R) = |e^'erfc(i?), (20) 
h{R) = ^[{l-2R^)IoR + ^R], 

Ira{R) - {-— + l-^)I,n-l{R) + —Ira-2{R) 

m Zm m 

where erfc(a;) is the complementary error function, 

erfc(a:) = j JTT^'^^' (^1) 

which rapidly decreases in x. {x)n — x{x -f l)(a:: + 
2)...[x + n ~ \) is the Pochhammer function. Several 
values of the complementary error function which can 
be obtained in any numerical library are: erfc(O) = 1, 
erfc(O.Ol) w 0.98872, erfc(O.l) w 0.88754, erfc(0.5) w 
0.47950, erfc(l)w 0.1573, erfc(2)w 0.00468, etc. 



A. Perturbative and exact calculations 

Following the e-I interaction whose matrix elements are 
expressed in Eq. (|18p . there is mixture between differ- 
ent states with different total radial quantum numbers. 
Consequently, the ground-state wave function e.g. for 
A'e = 1 electron is no longer the only one s- (n=0,Z=0) 
quantum state. Instead, it consists of several quantum 
states whose contributions Ci are different 

* = (22) 

i=l 

For Nf. = 2 electrons without impurity, by increasing 
A results in the ground state mixing of different config- 
urations. The two electrons will occupy higher Fock- 
Darwin states with decreasing energy spacing to lower 
their Coulomb repulsion. The easiest way to check the 
accuracy of the numerical results when the impurity is 
present is to consider the limitation when weak pertur- 
bation works perfectly. The first-order perturbation ap- 
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proximation: 



We 



E = Eo + ZJ2 < *o|l^(|r. - RDI^-o > +OiZ^) (23) 



where Eq is the GS energy of the QD with the ground- 
state wave function '^q in case without impurity is a good 
one when Z < 1. For the sake of simphcity in evaluating 
the perturbative part in Eq. for both = 1 and 

2 electrons, we assume that ^Pq in case Z = is a single 
term: the s-state. For Ne = 1 electron, this is always 
satisfied and its radial part is 



1 



(24) 
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FIG. 2: (Color online) Agreement between theoretically per- 
turbative (red dash-dotted) and exact diagonalization (black 
solid) calculations in the shift of the ground-state energy due 
to the presence of the impurity in single-electron (two upper 
curves) and two-electron (two lower curves) QDs for A = 0.1. 
The impurity effective charge is Z — —0.1. The inset shows 
addition data for the case that A is increased to 0.5 to exam- 
ine whether the perturbative calculations hold reliable. This 
energy shift is identical to the binding energy of the system 
and scaled to the Rydberg energy Ry = 5.31 meV for GaAs 
QDs. 

The situation changes for — 2. The above condition 
for to be a single term is only satisfied when A < 1. 
If so, the two electrons obeying Pauli exclusion principle 
with spins antiparallel will mainly stay in the s-orbital 
in case no impurity is present because their Coulomb 
repulsion is small as compared to the confining energy. 
The ground-state wave function is 

1*0 >= cl^cllO > . (25) 
Using this assumption, the total energy is estimated 



theoretically for the single-electron QD: 

-B(jVe=i) = ^0(„,=i, + 2^Z\nwoe^\Tic (i?) + 0(^2) 

(26) 

and the two-electron QD: 



4V^ZAfiwoe-"'crfc (i?) + O(Z^). 

(27) 

We plotted in Fig. [2] the energy shift due to the pres- 
ence of the impurity as a function of d using the above 
theoretical estimations (red dash-dotted) and the nu- 
merical results (black solid) for Z = 0.1 and A = 0.1 
{la = 1.03nm and fuvo = 200i?; = 1.062eV) for = 1 
(upper curves) and 2 (lower curves) electrons. Both the- 
oretical and numerical results are in good agreement. 

A small comment is made in case d is large. As seen 
from Fig. [21 the e-I interaction goes to zero slowly when 
d increases to a relatively large value, say d > a*g. The 



answer lies in the product 



' -erfc 



(«) 



where the ex- 



ponential function increases in d competitively with the 
complementary error function which decreases in d. 

The first-order perturbation theory works very well as 
long as Z < 1 and A is small enough. Because if A in- 
creases different e-I couplings of different configurations 
will occur which result in the presence of a substantial 
number of nonzero off-diagonal terms. This fact leads 
to an increasing difference between the first-order cal- 
culations and the numerical results. For example, in 
the two-electron QD as plotted in the inset of Fig. [2l 
this difference for i? = 0, which is also the largest 
value, is about 6% when A increases to 0.5. This 6% 
come from the configurations with minor contributions 
(ni = 1,^1 = 0;n2 ~ l.h — 0), (0,1;0,-1), etc. and 
n2 , 12 are the radial and azimuthal quantum numbers of 
the two electrons. 

Practically, one often has GaAs QDs with A > 1 i.e. 
larger QDs need to be considered. To describe such QDs, 
numerical calculations are used for different values of A. 
Here, we examine the case of A = 2 which has /q = 2a^ = 
20.6 nm and the corresponding hjjQ = R*/2 — 2.655 
meV. The Coulomb interaction unit (Vjf = 5.31 meV) in 
such a QD system is (2 times) larger than the confining 
energy hwQ. 



B. Singlet-triplet splitting energy 

The Coulomb interaction between the QD electrons 
strongly competes with the confining energy and with 
the e-I interaction as Z and/or A increase. Increasing A 
means that the confining energy becomes smaller with 
respect to the Coulomb interaction. As a result, elec- 
trons start to occupy higher Fock-Darwin levels. At zero 
B-field and in the absence of impurity, the ground state 
consists of several Fock-Darwin states where the dom- 
inant component is the s-wave term. In the presence 
of a charged impurity, those electron configurations that 
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FIG. 3: Splitting energy of the two lowest singlet and triplet 
energy levels as a function of (a) Coulomb interaction strength 
A for three different Z=-0.1 (blue dashed-dot), (black solid), 
and 0.1 (dash-dotted dot) and (b) impurity effective charge 
Z within the range (-1,1) for the case A = 2. The impurity is 
located at the origin of the single two-electron QD: d=0. 



fulfill L=const form different L-subgroups with nonzero 
contributions to the total wave function of the system. 
Increasing the effective charge Z results in a strong mix- 
ing between those subgroups. Consequently, there are 
more than 1 state which play as dominant components to 
the total wave function. We plot in Figs.[31[a) and (b), re- 
spectively, the splitting energy between the ground state 
(singlet) and the first excited state (triplet) as a func- 
tion of A and Z in case the impurity is located at the 
center of the single QD. As A increases [see Fig. [3l^a)], 
the splitting becomes smaller for both with and without 
impurity cases. It is because the singlet and the triplet 
states have many similar nonzero configurations. Details 
are found below from the discussion for particular val- 
ues of Z. From Fig. [3{a) we also notice that the energy 
splitting shift due to the presence of the impurity remains 
almost unchanged by changing A for both cases Z — —0.1 
(blue dashed dot) and Z = 0.1 (red dash-dotted-dot). It 
is because the impurity location is examined at the cen- 
ter of the QD d=0. Such splitting shift will change if 
the impurity is displaced to any other off-center position 
d^O. ' 

The Z-dependence of the splitting energy in Fig. ^h) 
for A = 2 shows a continuous decrease as Z changes its 
sign from negative (-1,0) to positive (0,1). The decrease is 
significant around Z=-l which coincides with the physics 
discussed above for the negative Z=-0.1 and -1 cases. 

It is found that for Z = —0.1 and —1 the ground state 
and the first-excited state in the two-electron QD con- 
taining a single charged impurity are the singlet and the 
triplet, respectively. Their major components are, re- 
spectively, the s-s and the s-p configuration, i.e. one 
electron is in the s- and the other in the s-(p-) orbital. 
Note that the s-s overlap has only one maximum at the 
origin. In opposite, the s-p overlap exhibits a minimum 
at the origin. 

With changing the position of the impurity, the energy 
spin splitting of the singlet-triplet states is obtained in 
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FIG. 4: Splitting of the two lowest singlet and triplet states 
as a function of the impurity position in a two-electron single 
QD for four different effective charges of the impurity (a) Z = 
-0.1, (b) -1, (c) 0.1, and (d) 1 in case A = 2 {hwo = 'i-K,!^ = 
2.655 meV). Dotted lines are the data obtained for the zero 
effective charge case {Z = 0) for reference. 



Fig- Ha) for Z=-0.1 and Fig. I^b) for Z=-l. The pres- 
ence of the impurity results in an increase in the splitting 
which is largest when the impurity is located at the cen- 
ter of the dot [as illustrated in Fig.lUJa) and Fig.HJJb) for 
Z = (dot lines) and for Z < (solid lines)]. 

The magnitude of the splitting over an impurity charge 
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/Z, with El the energy of 



the ground state and the first excited state in the doped 
QD, is found smaller for the larger Z case [see Figs.jH^a) 
and (b)] for a certain d value. This means that the e-I 
attraction becomes dominant over the e-e Coulomb inter- 
action. As a result, each state consists of different config- 
urations with compatible contributions to the total wave 
function. For example, by increasing the e-I effective 
strength from Z = -0.1 [Fig. Ufa)] to Z = -1 [Fig.H^b)] 
the contribution of the singlet s-configuration decreases 
from Cl « 0.83 to respectively « 0.64. The com- 
pensatory parts come from other configurations which 
also have L — such as (rti = 0, = 0; rt2 = 1, ^2 = 0), 
(0,0;2,0), (0,0;3,0), etc. Those states stay closer in energy 
with increasing A (= 2 in this case). 

Above we discussed the impurity effect for an attrac- 
tive impurity (positively charged Z <0). A negatively 
charged impurity which induces a repulsive coupling with 
electron is examined in Figs.|4l[c) and (d) for two effective 
charges Z=0.1 and 7j—1. The splitting energy between 
the ground state singlet and the first excited triplet is 
plotted as a function of the impurity position d. Ap- 
parently, when the impurity is found at the origin d=Q, 
it repels the two electrons most. The probability of the 
electrons to be found at the origin (s-orbital) reduces 
significantly. In this case, the electrons can be found 
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at other higher Fock-Darwin states. This means that 
the ground state energy becomes closer in energy to the 
first excited state. That explains a smaller splitting en- 
ergy we obtain for Z=0.1 [see Fig. SJc)] and Z=l [see 
Fig-IHd)] as compared to, respectively, the cases Z=-0.1 
[see Fig. ID^a)] and Z=-l [see Fig. ID^b)]. Moreover, we 
see from Figs. Hljc) and (d) that when the impurity is 
moved away from the center of the QD, the splitting en- 
ergy starts to increase. In other words, the probability 
of finding the electrons in the s-orbital increases. 



C. Impurity effect on the energy spectrum 

Technically, the presence of a charged impurity leads 
to an increase in the number of nonzero off-diagonal ele- 
ments of the Hamiltonian matrix. We present such exam- 
ination on the energy spectrum of the two-electron single 
QD as a function of the impurity position. The competi- 
tion between the two types of Coulomb interaction results 
in different relative orders of the energy levels depending 
on both the impurity charge Z and the impurity position 
d. 

For the weak perturbation case, Z=-0.1 [Fig.lSJJa)], the 
e-I interaction strength is much smaller (10 times) than 
the e-e interaction. The ground-state energy (level's de- 
generacy g=l) is the singlet state where the configuration 
with the two electrons in the s-orbital is highly dominant. 
The first- and second-excited states have degeneracies 
g — 2 due to the symmetry of L = ±1 states. Note 
that at B=0 the results are independent of spins. We 
present the results in the Sz=0 subspace. We discuss the 
L = 1 case. In the first-excited state, the most dominant 
configuration is (0,0;0,1), the second configuration is (0,- 
1;0,2), and the last configuration is (0,0;1,1). These three 
configurations and their exchange states have coefficients 
with opposite signs in the wave function (C(o,o:0,i) = 
— C'(o,i;o,o) ~ O-''', C'(o,_i;o,2) ~ — C'(o,2;o,-i) ~ 0-07, 
and the other C(i_i;o,o) = "^^(laiO.o) ~ 0.06) leading 
to the triplet first-excited state. The second-excited 
state has an opposite manner leading to the total spin 
5=0. Plus, the second dominant configuration and the 
third dominant configuration in the second-excited state 
switch their relative orders (sorted in probability) in 
the total wave function as compared to their orders in 
the first-excited state. The highest energies in Fig. [5] 
are the (L=±2,S'=0) states with the largest configura- 
tion (0, ±1; 0, ±1). Other considerable configurations are 
(0,0;0,±2) and (0,Tl;0,±3). 

In the two energy spectra presented in Fig. [Sj the 
ground-state and the first-excited state remain as singlet 
and triplet states. The fourth, fifth, and the sixth lev- 
els, which are the (L=0,S'=1) and (L=±1,S'=0) states, 
exhibit a crossing at d = O.OOOSa^. For d < O.OOOSa^j, 
the {L—0,S=1) state has a lower energy than the 2-fold 
degenerate (L=±l,5'=0) state. The major contribution 
to the wave function of the lower energy are, sorted in 
probability, (0,0; 1,0), (0, 0; 2, 0), and (0, 0; 3, 0). Increas- 
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FIG. 5: Low-level energy spectrum as a function of d of a two- 
electron QD in case (a) Z = —0.1 and (b) Z=-l for A = 2. 
Each level is labeled by: g as level's degeneracy, L, S as total 
angular momentum and total spin, and the major electron 
configuration (ni, Zi, 712, '2). and 712, '2 are the radial 

and azimuthal quantum numbers for electron 1 and 2. Inset 
in (a) is a magnification of the main plot which highlights the 
region with the occurrence of the crossing at d/a*g ~ 0.0008. 

ing d > O.OOOSa^ leads to the exchange between the 
(L=0,S'=1) and (L=±1,S'=0) states where the 2-fold de- 
generate (L=±1,S'=0) state has lower energy. 

The above physics is resulted from the dominant e-I in- 
teraction when Z=-l. The (L=0,S'=1) and (L=l,5'=0) 
states now have smaller energies than the {L=2,S=0) 
state. Besides, it is found that the 2-fold degenerate 
states (0, ±1, 0, T^), which dominate in the third-excited 
energy level {L=0,S=1) for Z = 0, have higher energies 
when the impurity is present. 

The positive-Z case does not affect the energy spec- 
trum as strongly as does for the negative-Z case. The 
reason is that the impurity in the positive-Z case induces 
the same type of Coulomb interaction with the QD e-e 
interaction. However, such an exchange in the relative 
order of the energy levels e.g. between the (L=0,S=1) 
and (2,0) states (corresponding to levels 4 — 6) can still 
be observed [see the crossing at d ~ O-Ma^ in Fig. [5]. 
The point d = 1.52a^ appears as a crossing-like point 
but it is only an almost-zero energy gap between the 
(1,0) and (0,1) states. Similar to the attractive impurity 
case, a larger e-I interaction also results in a smaller en- 
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FIG. 6: Low-level energy spectrum of two-electron single QD 
containing a repulsive impurity (a) Z=0.1 and (b) 1 for A = 
2. Inset in (a) is the magnification of the main plot into 
the region d=(0, l.la*g) to signify the crossing at d = 0.44a|3 
between the states (L=0,S=1) and (L=2,S=0). 



ergy gap between the first-excited state and the second- 
excited state. 



D. Summary for the impurity effect on the energy 
spectrum of 2-electron single QDs 

The cases of negative and positive Z affect the energy 
spectrum in different manners. Apparently, the impu- 
rity effect for the Z < case is expected to be stronger 
than the Z > case. When Z=-l the e-I interaction 
strongly competes and dominates over the e-e interac- 
tion. Let us take an example to illustrate the point. In 
the weak perturbation regime, i.e. \Z\ = 0.1, the e-e in- 
teraction dominates over the e-I interaction. The order 
of the low-level energy spectrum, then, exhibits no dif- 
ference between those two cases Z = ±0.1. Now, \Z\ is 
increased to 1. However, the system with negative Z ex- 
hibits a stronger e-I interaction. Such a strong effect for 
the attractive case Z=—l is observed in the exchange of 
the singlet and triplet states with L=l and L=0 as seen 
in the inset of Fig. [5jb). For the repulsive case Z—1, 
this exchange is no longer observed. Only the exchange 
between the triplet L=0 and singlet L—2 is found. The 
repulsive e-I interaction in case Z=l is manifested in the 



presence of the {L=2,S=0) state which was not seen in 
the lowest energy levels shown in Fig. \5lh) for Z=-l. 



IV. COUPLED QUANTUM DOTS 

Numerical results for the singlet-dot case were dis- 
cussed in detail mostly for A = 2 corresponding to 
Iq = 2a*g = 20.6 nm . This value of A was studies based 
on the realistic sizes of GaAs single QDs. However, we 
also theoretically discussed the results for a typical range 
of A=(l,2) as seen in Fig.[3]Ja). For the coupled QD prob- 
lem, we used the optimized parameter set after Ref.^ 
where = 10.01 nm {hujQ = 11.24 meV). Mapped on 
the single-dot case, the coupled 2-dot system will have 
"A"wl. The model of coupled dot system in the presence 
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FIG. 7: Schematic plot of the confining potential (bottom) 
and the charged impurity potential (top) in a coupled 2-dot 
system. We assume the impurity is negatively charged with 
Zi = 1 and is located along the line connecting the two 
double-well minima (±am) where ±, respectively, indicates 
the left (L)- and the right (R)-dot. The two solid lines on top 
depict the impurity potential on the two individual electrons 
in the two separate QDs. The upper solid line is the potential 
of the electron in the right- and the lower solid line of the 
electron on the left-dot, respectively. Dotted lines on top de- 
pict the other case when the impurity is located on the other 
half of the x-axis {xi < 0), i.e. it can be found inside the left- 
dot. Vb is taken to be 30 meV. The inset is the Coulombic 
potential of the impurity, located along the growth direction, 
which equally repels the two separate quantum dot electrons. 

of a charged impurity can be schematically described in 
Fig. [T] We plotted the impurity potential in the top 
panel for the case the impurity is located along the line 
connecting the two minima of the double well and in the 
inset for the case along the growth direction. It is clear 
that the impurity effect is largest when the impurity is 
at either of the two minima of the confining potential 
well in the former case. Whereas the latter case has the 
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largest impurity effect when the impurity is found at cen- 
ter of the system, i.e. R—0. We label the two impurity 
coordinates as Ri(2) = (a;i(2), 2/1(2), 2:1(2))- 

A. Singlet-triplet splitting 

1. Impurity position dependence 

In the coupled QD system without impurity and with- 
out magnetic field (5=0), the two lowest energy levels 
are the maximally entangled exchange spin states, re- 
spectively, the singlet ^tf and the triplet '^/^ . The next 
higher excited states are the linear combination of the 
two doubly-occupied singlets which result in "bonding" 
(^ ^2+^3 g^j^j "anti-bonding" ( '^^^^ ) states. 

Let us first consider the simplest case when the system 
contains only a single charged impurity Ri = (xi, yi, zi) 
(effective charge Zi) and the impurity plays only as a 
weak perturbation to the coupled dot system. However, 
there will be no restriction to the impurity location in or 
outside the system. Such a system allows us to provide a 
direct comparison to the single-dot case discussed earlier. 

In case the e-I interaction is attractive coupling, i.e. 
Zi < 0, the singlet-triplet exchange energy J is shown 
in Fig. EKa). When the impurity is located along the z- 
axis [see the black solid curve in Fig. [Hl[a)], the presence 
of the impurity increases the singlet-triplet spin split- 
ting J between the two electrons as compared to the 
case when no impurity is present [see horizontal grey dot 
curve in Fig. |Sl[a)]. This is understood as both electrons 
are attracted toward the impurity. Because the impurity 
equally couples to the electrons, the system favors the 
antiparallel electron spin state. This type of e-I coupling 
reduces the total energy of the system [negative binding 
energy presented as black solid in Fig. [HKc)]. As the im- 
purity is moved away from the origin [zi ^ - still the 
black solid curve in Fig.[5Ka)], J will decrease. Such a de- 
crease can be evaluated via Eq. (jA4p (see Appendix|X]) as 
the product of an exponential and complementary error 
function. 

When the impurity is located along the 'z'-direction 
of either the two separate dots, e.g. of the right-dot i.e. 
R-i = (am,0, zi) [see red dash curve in Fig. [Sfa)], J re- 
mains larger than the J for the case without impurity 
Zi — 0. However, J behaves very differently from the 
above two cases xi — and xi = Iq. The impurity 
not only no longer attracts equally the two QD electrons 
(similar to the xi = Iq case) but also affects the doubly- 
occupied states most. J in this case intersects the split- 
ting energy of the Zi = case at zi w Iq. Close to 
the QD sample, i.e. zi << Iq, J is always smaller than 
the splitting energy for the without impurity case, about 
10%. We reserve the detailed physical discussion around 
(±a,„ , 0) for a later discussion when we examine the case 
the impurity is located along the x-axis. 

We expect that a repulsive impurity induces an oppo- 



site spin order of the two electron spin orientations: the 
parallel spin state. This is illustrated in Fig. [SJb). A 
weakly repulsive impurity {Zi = 0.1), located along the 
z-axis (black solid curve), has the smallest exchange en- 
ergy (« 0.17 meV) when it is, apparently, at the origin: 
zi = 0. The reason is the impurity now repels both elec- 
trons. The e-I addition energy lifts up, few meV, the 
total energy of the coupled 2-dot system [see the posi- 
tive binding energy presented in Fig. [IJd)]. J rapidly 
increases and reaches the value of the non-impurity case 
as the impurity is engineered relatively far from the ori- 
gin, say 3^0 • 

The ground state consists of all three Hund-Mulliken 
singlets, i.e. the total wave function \E'gs — 
{^/if , V'f , V'f }■ However, two doubly occupied states play 
a small part, which are « 2%, to the total wave function. 

The binding energy of the impurity, which is defined 
as the energy difference of the system with and without 
a charged impurity, shown in Figs, and (d) was par- 
tially discussed above. In both cases Zi = tO.1, around 
zi = 0, the absolute value of the binding energy is found 
largest in case the impurity is placed closer to either of 
the two minima of the well (xi = a,„): 2.2 meV as com- 
pared to 2 meV and 1.53 meV of the xi — Iq and xi = 
case, respectively. Beyond a critical zi, say zi > 2lo 
the most dominant binding energy case (xi = am - black 
solid curve) becomes less dominant and compatible to the 
other cases xi = and ^o- AH three curves in Figs. [5{c) 
and (d) convert to the situation without impurity at the 
large zi limit. 

Now, we consider the case when the impurity is found 
inside the QD. In particular, we discuss the impurity ef- 
fect when the impurity is found on the x- or y-axis, along 
which the confining potential is constructed. 

Because the Zi = =f0.1 case was found to weakly affect 
the quantum dot qubits and that the doubly-occupied 
states have very small contributions to the singlet-triplet 
splitting, we can now use Heitler-London model to an- 
alytically check our numerical results. We found qual- 
itative agreement between the results presented in e.g. 
Fig. [5] and the analytical results shown in Eq. (|Bip . 
in particular the minimum (maximum) in J for Z=-0.l 
(0.1). Details are collected in Appendix IB] 

As signified in the top plot in Fig.[Tl the absolute value 
of the effective e-I coupling exhibits a maximum at the ei- 
ther of the two well minima and decreases rapidly as the 
impurity position is out of the minima [analytical ma- 
trix elements are presented in Eq. (jA5l) in Appendix |X] . 
However, the overlap between the two coupled QDs has 
a maximum at the origin (see e.g. Ref.^'^). In the weak 
impurity perturbation, i.e. Zi << 1, these two terms 
compete with each other. As a result, the singlet-triplet 
spin splitting J for Zi — —0.1 has the impurity posi- 
tion dependence as shown in Fig. O J has a maximum 
at the origin where the overlap between the two coupled 
dots is largest and the binding energy exhibits a maxi- 
mum around xi = l.6lo which is identical to the quasi- 
bottom positions of the left and the right dots ±a. As 
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FIG. 8: (Color online) Singlet-triplet splitting energy calculated as a function of zi coordinate of a single impurity in case Z\=- 
0.1 [(a) and (c)] and Z\ = 0.1 [(b) and (d)] for three different positions of the impurity a;-coordinate Xi=0 (black solid), xi^lo 
(red dash), and xi=am - either of the two double well minima (blue dashed-dot). VJ, is taken to be 30 meV. For comparison 
purpose, we recall the value of J ~ 0.204 meV as horizontal dashed-dot line in case no impurity is present for B—0. The insets 
(c) and (d) show the corresponding total binding energy of the three cases plotted in each main plot (a) and (b). 
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FIG. 9: (Color online) Singlet-triplet splitting energy as a function of the x-coordinate of the impurity position for (a) Zi — —0.1 
and (b) Zi —0.1 for 14=30 meV. The y- and ^-coordinates of the impurity position, Ry and Rz, are varied from the origin (black 
solid) to Ry = 0.2/o (green dotted), Ry = 0.5/o (violet short-dashed), Rz — 0.2lo (blue dash-dotted), Ry = 0.2lo,Rz ~ 0.2Zo 
(red dashed), and R^ — 0.5Zo (dash-dot-dotted). The omitted coordinates are implied to be 0. We cover many different possible 
positions of the impurity such that: it can be along the a;-axis or along the y-axis, or around the either the two minima with 
both xi and yi coordinates slightly changed while zi — 0, etc. Notice the case where the impurity is shifted along the y-axis. 
This shift does not add much physics to the system. As for example, see the insubstantial difference between the exchange 
energies in the black solid (yi = a;i = 0) and green dot (j/i = 0.2io) curves, or between the blue dashed-dot (z\ = 0.2Zo) and 
red dash (j/i = z\= 0.2io) curves, etc, presented in both (a) and (b). 




compared to the case the impurity is located along the 
^-direction, the considered case has a minimum around 
the point x\ « 2Zo (see Fig. [3]) which is identical to the 
analytical minimum obtained using Heitler-London ap- 
proximation (see Appendix |B|. We note that such a 
minimum is obtained for an attractive impurity which 
has a positive charge Zi < and \Zi\ « 1. Moving 



the impurity out of the dot system (blue dashed-dot, red 
dash, and cyan dash-dotted dot in Fig. ^ will lead to a 
decrease in J. This situation can be considered as the 
coupled QDs interacting with charged impurities found 
close to the surface during the growth process. In partic- 
ular, the case that the impurity y-coordinate is displaced 
to < yo/^o < < 1 has only small difference in J from the 
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case the impurity is positioned exactly along the x-axis, 
Uo = (see the green short-dot curves and solid black 
curves in Fig. [9]). 
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FIG. 10: (Color online) (a) Singlet-triplet spin splitting J 
as a function of the impurity effective charge Zi within the 
range Zi = ( — 1, 1) for three different positions of the impurity 
along the z-axis and (b) the singlet-triplet state-diagram as a 
function of the impurity charge Zi and its position along the 
z-axis. Vi,=30 meV. S and T refer, respectively, to the singlet 
and triplet state. When the impurity is positioned close to 
the origin [zi^O (black solid) and zi=0.25Zo (red dash) in (a)], 
the system exhibits a triplet-singlet transition at Zi — 0.54 
and 0.8, respectively. Green dot curve is added for reference 
purpose which shows that as the impurity is found at (or 
close) to either left- or right-minimum of the well the system 
is always the singlet. 



2. Impunty charge dependence 

The above critical point xi ~ 2lo at which the 
singlet-triplet splitting exhibits a maximum/minimum 
(see Fig. ^ depends on the effective impurity charge Zi 
and the well barrier height Vb- In this subsection, we ex- 
amine the Zi dependence of the singlet-triplet spin split- 
ting J. 

Besides, we notice the physics sampled around the ori- 
gin in Fig.[8jd) and [DJb) where the presence of a repulsive 
charged impurity significantly lowers the singlet-triplet 
spin splitting J. It is expected that by further increasing 
the impurity charge Zi, the system can be visited in the 
triplet state, i.e. J < 0. 

Such triplet-singlet transition occurs as Zi is larger 
than w 0.55 in both cases the impurity is located along 
the z-axis and along the x-axis as shown, respectively, 
in Figs. [To] and [TT] Our calculations for the attractive 
coupling case Zi < show that there is no triplet-singlet 
transition (see Figs.lTUland [TT|) due to the fact that a pos- 
itively charged impurity attracts both electrons therefore 
the favored spin state is always the singlet. 

The triplet-singlet transition is further explored for dif- 
ferent negative charges Zi > as seen in both Figs.fTUTa) 
and lllf a) for different impurity locations: at the center 



[black solid in Fig. [TOfa) and full rectangles with black 
solid in EJa)] and off-center - slightly away from the 
origin. We only examine the triplet-singlet transition for 
the systems with the e-I coupling smaller or compatible to 
the e-e interaction. In fact, the limit of Zi = ±1 does not 
bear much physical meaning. However, we theoretically 
examine that limit to support a complete understand- 
ing of the effective e-I strength on the exchange electron 
qubits in zero B-field. 

It is undoubted that when the impurity equally in- 
teracts with the two quantum-dot electrons and the e- 
I interaction strength is considerable the triplet ground 
state occurs at a smaller Zi , as compared to the unequal 
electron(s)-impurity interaction case. This argument is 
clarified in e.g. Fig. [TIT a) where the singlet state occurs 
by increasing Zi to w 0.54 for i?i = and further up to 
0.7 for i?i = a,„/4 (w 0.43^0) • On the other hand, at a 
certain strong enough Zi , the triplet state in case the im- 
purity is located along the x-axis remains longer than in 
case the impurity is located along the z-axis. For exam- 
ple, for Zi = 1, we obtain that the triplet state stays up 
to xi sa 0.6Zo while it is found only up to zi « 0.46^o- It 
is worth noting that at a higher Vb = 35 meV, due to the 
influence of the impurity potential, the triplet state re- 
mains up to a larger xi in comparison to the 30 meV case. 
For example, such a triplet-singlet transition is obtained 
at xi — 0.42^0 and 0.72^o for Zi—0.6 and 1, respectively. 

When the impurity is located at the right well bottom 
for Zi — =pl, the ground-state is always a singlet however 
the maximally entangled component tpf is replaced by 
one of the two double-occupied components {V'f j'^f } as 
the major contribution to the total wave function. 

We summarize in Figs. fTUT b) and [TlTb) the occurrence 
of the triplet or singlet as the ground state when changing 
the impurity position i?i and its effective charge Zi for 
Vb = 30 meV. The triplet state starts to occur when 
Zi is increased to « 0.54 at Ri = 0. The largest triplet 
region = (0, 0.38/o) in Fig.[TOi;b) and (0, O.G^o) in Fig. 
[TTTb)] is, apparently, seen for the largest considered Zi 
{—1). In particular, close to the origin, the triplet region 
is rapidly shortened. For example, for Zi=0.6 the triplet 
region is in Ri — (0, 0.3) while such a region squeezes to 
(0, 0.07) as Zi decreases to 0.54. Note that we show the 
singlet-triplet state-diagram on the right half of the x- 
(2:-)axis. As the impurity is positioned along the other 
half of the x- (z-)axis, the state-diagram is found similar. 



3. Inter-dot separation dependence 

In this subsection, we study the inter-dot dependence 
of the singlet-triplet splitting J. First, we place the im- 
purity at the center of the dot system, i.e. Ri = (0, 0, 0), 
and tune the barrier height from Vb — 35 meV down 
to 13 meV. In the meantime, the inter-dot separation 
(2am) will decrease from 2am = 35 nm (« 3.5^o) to 17 
nm (ss 1.6/0) ■ It is worth noting that the case that the 
impurity is placed exactly at either of the two well bot- 
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FIG. 11: (Color online) (a) Singlet-triplet spin splitting of a coupled 2-dot system as a function of the effective charge Zi within 
the range —1 < Zi < 1 for xi = (full squares with black solid) and xi = am/4: (full circles with red solid), (b) Singlet-triplet 
state-diagram plotted in the impurity-position— effective-charge _Ri — Zi plane. Vt is taken to be 30 meV. Horizontal dot line in 

(a) is used to clarify the triplet-singlet transition. In the region close to the center of the system, at ~ 0.07/o (about 0.04am), a 
relatively large repulsive impurity potential, say Zi = 0.54, is enough to induce a triplet-singlet transition [see the state-diagram 

(b) ]. The stronger the effective charge, the broader the triplet region. Consequently, when moving the impurity close to the 
right bottom it requires a larger Zi to observe the triplet-singlet transition occur. S and T stands, respectively, for the singlet 
and triplet state. 




FIG. 12: (Color online) Singlet-triplet splitting energy J cal- 
culated as a function of the inter-dot separation 2am for dif- 
ferent positions of the impurity when it is located along the 
2;-axis for Zi=-0.1, Z2=0. We added data in case no impu- 
rity is present as the open squares with magenta dash line. 
Because 14 changes (from 35 meV down to 13 meV) effec- 
tively controlling the inter-dot separation, i.e. am, while the 
impurity position should be fixed during the Vt, modification, 
we examine different locations of the impurity. Several ex- 
amples of corresponding (Vf,,2a„i) are (35 meV,3.5/o), (30 
meV, 3.3Zo), and (25 meV,2Zo). Inset is the magnification of 
the region sampled by the big open circle for high-Vi,-limit. 



toms does not have any physical meaning in this case 
because the bottoms of the double-well change concomi- 
tantly with Vb changing (see Fig. [Ij. In fact, the im- 



FIG. 13: (Color online) The same plot as Fig. [12] for Zi = 
0.1. As the impurity position is engineered further away from 
the origin, the double-dot system reveals crossings between 
different J for different _Ri. 



purity position should be fixed thus its position only 
varies relatively with respect to the double-well minima 
of different- Vf, systems. Consequently, different- Vf, 2-dot 
systems exhibit various impurity effect on J. We obtain 
the inter-dot separation dependence of the singlet-triplet 
spin splitting J for Zi = —0.1 in Fig. [T2] and Zi=0.1 in 
Fig.m 

For Zi = —0.1, the splitting energy J becomes smaller 
as the inter-dot separation increases (see Fig. [T^ . How- 
ever, J is found larger than that of the case without 
charged impurity as illustrated in e.g. the full rectan- 
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gles with black solid line (i?i = 0) and open rectangles 
with magenta short dash (no impurity) in Fig. 1121 In 
this case, this means that the potential well height is 
"weakened" by the attractive impurity which attracts 
the electrons toward and therefore supports the anti- 
parallel spin interaction of the two electrons. Displac- 
ing the impurity away from the center of the double- 
dot system (full rectangles) along the x-axis results in 
changes as seen in a series of J presented as full cir- 
cles with red sohd (for i?i = 0.5^o), full up-triangles 
with green dot (for i?i = 0.7/o), full down-triangles with 
blue solid (for i?i — Iq), full rhombuses with dashed-dot 
(for Ri — 1.5Zo), full right-triangles with yellow short- 
dot (for Ri = 2lo), full circles with dash-dotted-dot (for 
Ri = 2.5Zo), and full stars with solid (for Ri — 31q). 
First, we analyze the data for the lowest case of VJ, = 13 
meV (corresponding to 2am ~ 1.65/o)- The minimum in 
J for this case is found at xi « while the maximal 

value in J is always found at the origin i?i = 0. Moreover, 
the system tends to covert to the no-impurity situation 
when the charged impurity is engineered far enough from 
the center of the dot system. Therefore, from the i? = 0- 
case to the R = 0.7/o-case we obtained such a decrease 
in J from « 2.12 meV to 1.87 meV but for R = Iq we 
obtained J almost identical to the J of a similar system 
but without impurity. Note that this impurity position 
Iq is considered still close to the right-bottom of the well 
barrier. This feature is similar to the physics around the 
impurity position R w 1.2/o for Vt = 30 meV as previ- 
ously examined in full squares with black solid in Fig. [HI 
Further increasing xi to 1.5lo and 21q (as seem in the full 
rhombuses with cyan dashed-dot and full right-triangles 
with short-dot, respectively) we obtain the increase back 
in the latter case. This is found due to the fact that both 
1.5/o and 2lo are located on the right-hand-side of the 
minimal point in J. The latter case exhibits an increase 
because the system to this extent tends to convert to 
the non-impurity case. In the furthermost case Ri = SIq 
(full stars with solid line) J is found closer to that of 
the non-impurity case which has higher energy than the 
cases Ri = l,5lo,2lo, and 2.5Zo. 

Such a shift in the maximum of J by changing Vf, can 
be seen in Fig.[T2]and its inset for the high-Vf, limit where 
the furthermost impurity positions Ri = 1.5/o, 2/o, 2.5Zo, 
and 3^0 stay closer to the non-impurity case in compar- 
ison to the small- Vf, cases. This is related to the fact, 
which was mentioned above, that the relative distance of 
the impurity position to the minimal point of J varies 
where the impurity can be found either in the left- or 
the right-hand-side of the minimal J as Vt changes. For 
Vb — 30 meV (corresponding to 2am ~ 3.3/o), from 
the center-i?i — 0- to the /o-case the singlet-triplet 
spin splitting stays higher in energy while the others 
(i?i = 1.5lo,2lo,2.5lo, and 3/o) have lower J than the 
non-impurity case. The minimal point for this system 
is 2/o- Therefore, the i?i = 1.5Zo-case has a higher J 
than the R = 2Zo-case which was opposite in the previous 
case Vb = 13 meV. As a consequence, there appear cross- 



ings in J for different Obvious crossings are seen 
for the cases Ri = 1.5lo,2lo,2.5lo, and 3lo- For exam- 
ple, J for the Ri = 0.5Zo-case (full rhombuses with cyan 
dashed-dot) stays lower in energy however it changes as 
Vb varies and is found higher (or equal) in energy than 
the J in the others (see the inset as a magnification for 
the high-Vfc-limit). 

Apparently, when placing the impurity close to the ori- 
gin we always obtain a larger singlet-triplet spin splitting 
as compared to the non-impurity case. The largest case 
has AJ « 0.6 meV for i?i = and V& = 13 meV. It is 
worth noting that a larger Vb results in a smaller differ- 
ence in J which can be understood as a lower "tunneling" 
rate of the two electrons in the two separate QDs. Quan- 
titatively, such a decrease can be roughly estimated by 
Eq. (|BT]) 

It was made clear in the presence of an attractive 
charged impurity the e-I coupling increases the "tunnel- 
ing" rate between the two dots. Positioning the impurity 
at different locations inside the coupled 2-dot system, 
in particular around the {xi « 2/o,0) point, gives rise 
to different Vt,-dependent singlet-triplet splitting which 
has crossings between different i?i-curves. As illustrated 
now in Fig. [T2]for Zi = 0.1 case, J in the i?i = case 
(full squares with black solid) stays lowest in energy and 
the J-difference as compared to the non-impurity case is 
found largest. For example, AJ « 5.5 meV for Vf, — 13 
meV. J in the furthermost case i?i — 3lo (full stars with 
navy solid) stays slightly above the non-impurity curve 
(open squares with magenta short-dot). This implies that 
such an impurity position is on the right-hand-side of the 
maximum point in J for all the considered Vb [g(13, 35) 
meV] in Fig. [T^l It is found opposite to the Zi — —0.1 
case where the i?i — 3lo-J-cuive stays slightly lower 
than the non-impurity one. Crossings between different 
i?i-dependent J-curves are still obtained in the current 
case Zi = 0.1, e.g. the one around 2aTO = 2.1/o (corre- 
sponding to Vb = 15.5 meV) between the Iq- (full down- 
triangles with blue solid) and the 3Zo- J-curves, or the one 
at 2am ~ 2.6I/0 (Vb = 20 meV) between the 1.5/o- (full 
rhombuses with cyan dashed-dot) and the 2/o-J-curve 
(see also the inset of Fig. [T5)) . etc. 

It is important to note that the \Zi\ —0.1 case studied 
so far has not exhibited any crossing between the sin- 
glet and triplet states, i.e. the triplet state always stays 
higher in energy than the singlet state. The presence of 
an attractive impurity increases the singlet-triplet split- 
ting opposed to the effect seen for a repulsive impurity. 
The extreme behavior of J is obtained at either of the 
two minima of the potential well. An attractive impu- 
rity seems to increase the "tunnel" rate of the coupled 
double-dot system. 



B. Energy spectrum 

In the single QDs as discussed in Sec. Ill, impurity 
effect on the energy spectrum of the system is substantial 
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FIG. 14: (Color online) Energy spectrum of the QD system 
studied (a) in Fig. [SJa) and (b) in Fig. |5^a), respectively, as 
a function of the impurity position (along the z and the x- 
axis) for Z\ — —0.1, Z2 — 0. Vt = 30 meV. Because there is 
no crossing/anti-crossing in the energy spectrum, we use the 
same style of line for a certain 7?i-dependence. In (a), solid 
curves depict the impurity-position dependence when the im- 
purity is located along the z-axis and dash-dotted curves de- 
pict the situation when the impurity is outside the double- 
dot system [along the line {yi = 0; a;i = 0.5/o)]. In (b), solid 
curves depict the case _Ri is changed along the a;-axis and 
dash-dotted curves depict the situation where the impurity is 
outside the double-dot system but along the i-direction. The 
latter case of (a) with xi = 0.5lo (dashed-dot curves) reveals a 
stronger impurity effect on the coupled qubits in comparison 
to the case when the impurity is on the z-axis [solid curves 
in (a)]. The Ri = (0, 0, z;i)-case induces equal exchange cou- 
pling to the two electrons mostly found on the individual dots. 
The latter case of (b), in opposite, has much weakened im- 
purity effect on the double-dot system in comparison to the 
case where the impurity is found on the line connecting the 
two confining potential minima. 



as e-I interaction is competitive to the e-e interaction. 
We obtained different crossings betvireen lovif-lying excited 
states and anti-crossings with varying energy gaps. If 
such a similar effect is found in the coupled double QD 
system, the question 'Vhether quantum operations in the 
coupled dot system are affected" attracts our attention. 

In order to examine the validity of using coupled dots 
containing charged impurities as the basis of quantum 
logic gates in quantum computation, we mainly look into 
the low-level energy spectrum. 

The energy spectrum of double QDs containing a sin- 
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FIG. 15: (Color online) The same plot as shown in Fig. [14] 
for a repulsive impurity Zi — 0.1. There is a change in the 
relative order of the two excited states in (b) as compared to 
the above case Zi = —0.1 (see Fig. I14[). 



gle impurity located at Ri = (0,0, zi) and (0.5Zo, 0, zi), 
respectively, is shown as solid and dash-dotted curve in 
Fig. [Hf a). The latter case (dash-dotted) has a larger 
impurity effect on the energy gap between the first two 
lowest levels and the two highest energy levels which can 
be understood in an earlier discussion. The energy levels 
bend down to the zi = as its s-coordinate xi ^ as 
compared to the case a;i = because the impurity has a 
stronger interaction to the electrons when it is closer to 
the well bottoms. 

The xi-dependence of the energy spectrum is studied 
in Fig. [TW b) for two different z\ — Q and z\ — 
The minimum, discussed earlier in Fig. [S] related to the 
physics of having an impurity around the well bottoms, 
is seen in the energy levels of Figs. ITW b) for zi — 0. The 
spin splitting between the two fully-filled states is most 
affected also when the impurity is at the well minima. 
The probability of finding the two electrons in the right 
dot (containing the impurity) increases. The absolute 
value of the e-I coupling for the right (in the present 
case) doubly-occupied state is maximal at (am,0). 

For Zi — 0.1 [see Fig. [T5l b)]. it is found opposite to 
the attractive Zi = —0.1 case. The impurity (in the 
right-dot) now tends to repel both electrons which leads 
to a higher probability of finding the two electrons in the 
left-dot. Consequently, the doubly-occupied state of the 
right-dot gains energy via the e-I interaction and stays 
higher in energy (blue dashed-dot) than the singlet state 
(green dot) of the left-dot. Such a maximum around the 
bottom {am , 0) in the blue dashed-dot curve in Fig. [T5l b) 
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has similar physics discussed earlier for Z\ = —0.1. Re- 
meraber that in the attractive case, the singlet-triplet 
spin splitting has a minimum around 2/o (see Fig. [9]) 
which is resulted due to i) the competition between the 
wave function overlap of the coupled 2-dots and the impu- 
rity potential on the right-dot electron and ii) the doubly- 
occupied states have relatively small contribution to the 
total singlet state (the ground state). Now, the effective 
charge changes its sign (0.1) and we obtain a maximum 
[see the inset of Fig. lTST b)]. The other left-dot fully occu- 
pied state [green dot in Fig. fTSl b)] only shows a decrease 
in x\. 

Even though the presence of the impurity changes the 
singlet-triplet splitting as discussed above, the energy gap 
Ae between the highest spin state that stores information 
and the lowest unwanted state, in this case the first and 
the second excited states, remains much larger than the 
singlet-triplet splitting energy J . In case Vf, = 30 meV, 
J/Ae is typically < 0.07. This means that adiabatic 
condition is satisfied so that higher excited states are not 
involved when the system is evolved to its desired state. 



C. Strong electron-impurity coupling destroys 
coupled qubits 

1. Strong perturbation 

In the following, we discuss in detail the influence of 
a strongly perturbative impurity on the energy spectrum 
for Zi = =f1. Results for other intermediate Zi (e.g. 
=F0.6, =f0.8) are collected in Appendix ICl 

We plot in Fig. [1^] the impurity-position dependence of 
the singlet-triplet spin splitting J when the impurity is 
located along the z-axis for (a) Zi = — 1 and (b) Zi — 1. 
The former case has a similarly qualitative curve to the 
case Zi = 0.1 with a much larger exchange energy J. The 
latter case has a triplet-singlet transition which occurs 
at zi = 0.38^0 as obtained already in the singlet-triplet 
state-diagram Fig. [TUfb). This property is resulted from 
the strong e-I interaction as Zi is increased to 1. 

For the case the impurity is located along the x-axis 
and Zi — 1, we obtain in Fig. 1171 and its inset the entire 
triplet state for the case the impurity is engineered at the 
origin or very close to the origin, say i?i < 0.6lo (see the 
inset). The i?i = 0-case has the impurity that equally 
repels the two electrons in the two individual QDs and 
the impurity is kept distant from the two well minima. 
As a result, the favored state becomes the triplet with 
two electron spins aligned parallel to each other. Moving 
the impurity off-center means that there appears a bias 
in the impurity-electron coupling with the two electrons. 
Such a triplet state becomes weakened and disappeared 
when the bias increases to its maximum for the case the 
impurity is found around the bottoms of the confining 
potential. 
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FIG. 16; (Color online) Singlet-triplet spin splitting J calcu- 
lated as a function of the impurity position along the z-axis 
for (a) and 1 (b). Vt = 30 meV. S and T stand for, 

respectively, singlet and triplet. The horizontal dotted line 
in (b) is used to clarify the triplet and singlet states. Such a 
triplet-singlet transition is obtained at zi = O.Slo for Zi=l. 




FIG. 17: (Color online) Singlet-triplet spin splitting J of a 
singly-doped charged impurity coupled double-dot system as 
a function of the impurity position when it is located along the 
a;-axis for Zi = 1. Vi, = 30 meV. Inset is the magnification 
of the circled region in the main plot which highlights the 
triplet-singlet transition in the small xi region. 



2. Destroy of coupled qubits 

Impurity effect of small effective charges shown in 
Figs. [9] and [15] as discussed above for Zi = =f0.1 results 
in a relatively small coupling between the three Hund- 
Mulliken singlets in ground-state singlet. Quantitatively, 
the major (w 98%) component is the V'f therefore the 
mixing of '0f with the other two singlets "^f 3i with the 
probability about 1% each, can be neglected. The consid- 
erable impurity effect is observed in the excited singlets 
where such full bonding- and anti-bonding states are no 
longer found at B=0 when the impurity is found in either 
of the two individual QDs. Plus, the probability of find- 
ing the electrons in a fully-occupied state is increased. 

In the next examination, we study the case when the 
coupling between the electron and the impurity is com- 
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FIG. 18: (Color online) Energy spectrum of a coupled QD system exhibits strong impurity effect examined as a function of the 
(single) impurity position located along the a;-axis for a (a) positively {Z1—-I) and (b) negatively (Zi — l) charged impurity. The 
inset in each layer is the corresponding singlet-triplet spin splitting J where both exhibit a maximum around the right-bottom 
of the well. Different from the attractive case Z\ = —1 m (a), the repulsive case Z\ = 1 m (b) overcomes a triplet-singlet 
transition at « Q.Qlo [see the inset of (b)]. The green dot curve in (a) represents either the -i/jf singlet or the tpf singlet (see 
discussion in text) while the green dot in (b) represents either the ipi singlet or the -^f singlet. Red dash in both (a) and 
(b) is the triplet state. Circles that sample the anti-crossings indicate the transition of i/)f as the ground-state to the first 
excited state and the transition back to the ground state as x\ increases. Anti-crossings and almost zero energy gap between 
the first and second excited states in both (a) and (b) are the evidences to claim that the coupled qubits are "destroyed" and 
the adiabatical condition is violated. 
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FIG. 19: (Color online) Probability of finding the three singlet 
states as defined in Eq. (I17p in the ground state as a function 
of the impurity position for different Zi for weak [Zi=-0.1(a) 
and 0.1(b)] and strong [Z\=-l{c) and 1(d)] impurity potential. 
The triplet state can be the ground state and only "couples" 
(e.g. via an anti-crossing as discussed earlier in Fig. I18|l but 
does not mix with the singlet states therefore the probability 
of finding such state is always 1. The impurity position is 
Ri = (a;i, 0, 0). Top dash-dotted horizontal line corresponds 
to the probability of finding the state in case without 
impurity for B=0. 



patible to the electron pair Coulomb interaction by in- 
creasing the absolute effective impurity charge Zi to ±1. 
The energy spectrum are now shown in Figs. llSf a) and 
(b). For the attractive impurity case [see Fig. [TSla)]. 
the triplet remains as the first excited state while the 
ground state now mixes three singlets with the leading 
term varies either the spin-exchange singlet, , or the 
right-dot doubly-occupied state f/jf . The other singlet 
-01 (left-dot) has slightly larger than zero contribution. 

Particularly, the singlet-triplet splitting in case Zi = 
— 1 changes mostly differently in comparison to the weak 
impurity potential case Zi — —0.1. J now has somewhat 
similar physics to the case Z = \ around the right bottom 
of the potential well. A maximum in J (« 6 meV) was 
seen on both cases as shown in the insets of Figs. [TSFa) 
and (b). 

Now, the strong e-I interaction leads to the triplet- 
singlet transition at i?i ~ O.6/0 (as discussed already) for 
Zi = 1 and "virtual" coupling between the three states 
"0"^, V-'f J and 02 with each other. 

As compared to the weak perturbative case where the 
minor contributions of the ip2 3 states, up to 1% each, 
are negligible, such contributions now start increasing 
and can exceed 50% and even larger for the doubly occu- 
pied singlet of the dot containing the attractive impurity 
(the right-dot i.e. the ip§) and of the dot without im- 
purity (the left-dot i.e. the -^f ) for the repulsive case. 
Consequently, the doubly-occupied state can become the 
ground state and the exchange spin singlet becomes the 
first excited singlet. From Fig. rTST a). the green short 
dot corresponds to the (major) and the blue dashed- 
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dot the ^/if (major) state in the region Ri < O.SIq and 
Ri > 2.7lo. The region in between O.8/0 < Ri < 2.7Zo 
has the ground state as the doubly occupied singlet ^/if in 
Fig. [TST a') and i/'f in Fig. [TST b). The minor components 
in the first excited state in both cases, respectively, V'f 
plus and V'f plus V'f i fluctuates around 7% and 3% 
as the impurity is located close to the center of the sys- 
tem. In Fig. [TSlb') we see a switch between the left- and 
right-doubly-occupied singlets where the singlet of the 
right-dot (with impurity - blue dashed-dot) stays higher 
in energy than the one of the left-dot (green dot). 

The "virtual" coupling between the three singlets and 
the triplet addressed above is now examined in both cases 
Zi = =f1. At the well right-bottom i?i — (a,„,0), the 
triplet (red dash) "couples" with the first excited sin- 
glet as green dot curve [V'f in Fig. rTST a)] and [-^f in 



Fig. [Wb)]. The energy gap between the triplet and the 
first excited state is found close to zero. This virtual 
coupling can be understood by the two anti-crossings be- 
tween the tpf and in Fig. fTSla) and between tpf and 
V'f in Fig.IIHJb). 

To better understand the physics manifested in the 
mixing of different singlet states as well as the "cou- 
pling" between the triplet and the first excited singlet 
as shown in Figs. ITSTa) and (b) and the triplet-singlet 
transition in case Zi = 1 we calculate the probability of 
finding the three singlet states as a function of the impu- 
rity position inside the double-dot system for weak and 
strong e-I coupling, respectively, in Figs. [T9l a) and (b) 
and Figs. \T9\ c) and (d) for the ground state. We first 
explain for the triplet-singlet transition in case Zi = 1 
[see Fig. [TWa)] and the "coupling" of the first excited 
singlet and the triplet in case Zi = =f1 [see Figs. [T8l a) 
and (b)]. The triplet-singlet transition can be under- 
stood from Fig. [TW d) where the region xi/Iq < 0.6 has 
zero probability of finding the V'f state in the ground 
state. It is because the triplet state becomes the ground 
state and the singlet ipf is the major component of the 
first excited state wave function in this considered region. 
The dominant of the doubly occupied singlet V'f in the 
ground state for Zi — ~1 and of the V'f for Zi = 1 is now 
specified as the region O.SZq < xi < 2.7/o in Figs. ITWc) 
and (d). This region is identical to the maximal behavior 
of J as seen in Figs. [TST a) and (b) and their two insets. 
For weak impurity potential cases Zi ~ ^0.1 there are 
no anti-crossings in the energy spectrum and the doubly- 
occupied singlets remain always as the excited singlets 
with small probability [see Figs. iroT a) and (b)]. 



D. Coupled qubits perturbated by two impurities 

For simplicity, we first assume that two impurities have 
identical charges Zi = Zi and one impurity is kept at one 
of the two minima R2 = (— am,0, 0). The other impu- 
rity can be arbitrarily located along the other side of the 
X-axis, Ri = (xi > 0,0,0). Under these circumstances, 
the two impurities generally induce different e-I couplings 
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FIG. 20: (Color online) Singlet-triplet splitting (black solid) 
of a coupled two-dot system with the presence of two identical 
charged impurities but located asymmetry along the i-axis: 
one is fixed at the left bottom and the other is located along 
the X > 0-axis for weak e-I coupling (a) Zi = = —0.1 and 
(b) Z\ = Z2 = 0.1. Vh is taken to be 30 meV. Violet dashed- 
dot curves are the J for the single-impurity system extracted 
from Fig. 



and hence affect the qubits asymmetrically except when 
xi = — X2 = a,„. In this case, both doubly-occupied 
singlets are found compatible which have similar contri- 
bution to the two excited singlet states. Therefore, we 
no longer see a minimum or maximum, depending on the 
sign of Z\ , in the second or third excited state as obtained 
in case Zi — 0. Such a pronounced minimum is now 
"eliminated" due to the presence of the second impurity, 
located at the other bottom (left) of the 2-dot system, 
which induces a competing e-I coupling. The singlet- 
triplet spin splitting J has maximal behavior when the 
impurity is found close and at the bottom of the right- 
dot. 

Another way to study the asymmetric impurity effect 
in the considered double-dot system with two impurities 
is to adjust the effective impurity charges with respect to 
each other while their positions are kept symmetrical to 
the origin, i.e. ZijZi ^ 1 while Ri=-R2. We compute 
the singlet-triplet spin splitting for three different values 
of Zi/Zi = 1,1.1, and 1.2 while Ri^i changes for both 
positively and negatively charged cases in Fig. [2lT a) and 
(b), respectively. The results shown in Fig. [^l] are ob- 
tained for Zi/Zi slightly different than unit to examine 
the asymmetry arisen due to the dissimilarity in impurity 
charges as compared to the asymmetry due to the impu- 
rity positions discussed earlier in Fig. [501 for Zi — Zi. 

Apparently, if one replaces one of the two identical im- 
purities by an opposite charge to the other one can revisit 
the "asymmetric" phenomenon between the two doubly 
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FIG. 21: (Color online) Singlet-triplet spin splitting J as a 
function of the impurity position in case the two impurities are 
symmetrically located, Ri=-R2, but they are charged slightly 
differently (a) both impurities are positively charged and (b) 
both impurities are negatively charged with Zi = ±0.1 and 
Z2 = Zi (black solid), l.l^i (red dash), and 1.2Zi (blue dash- 
dotted). Vb is taken to be 30 meV. Vt is fixed at 30 meV which 
corresponds to am ~ 1.64/o. Note that the system converts to 
the situation without impurity at a much larger xi than the 
singly-doped impurity case (previously seen around 3.5Zo)- In 
this plot, at xi — Alo the system does not convert to the case 
without impurity (e.g. A J ~0.01 meV) and the total binding 
energy AE ^0.8 meV for Zi = -0.1 [black solid in (a)]. 



occupied singlets in the two dots which now can be en- 
hanced as a double. Our further calculations show that 
if Zi=-Z2=-0.1 and Ri = -R2 = (-0^,0,0), for V&=30 
meV the two doubly-occupied states have a significantly 
large energy gap and the higher level, which is the right 
occupied state in this case, exhibits a maximum around 
(omjO) - one of the two potential well minima. Gen- 
erally, this considered situation can be extended to the 
case that the charged impurities are different types and 
obey a specified charge distribution. By examining the 
energy spectrum of the coupled dot, as illustrated in the 
above discussion, one can detect the impurity effect on 
the singlets of the double-dot system. This can be made 
possible due to the distinguishability between the total 
positive and the total negative charge cases. It is also 
possible that both impurities are found either in the left- 
or right-dot, however, not much different physics from 
the case with only a single impurity will be added. 

The last part of this section is reserved to examine 
one special case described in what follows. The cou- 
pled 2-dot system interacts with two negatively charged 



impurities Z\ = Zi = 1, equal to the electron charge. 
Now, the two impurities are engineered such that they 
are far enough from the target (A) coupled 2-dot sam- 
ple and only induces electrostatic interaction to the two 
electrons of the target system. We examine the case 
that the two impurities have varied positions but their 
center-of-mass coordinates are kept unchanged, in this 
case R^2 = (Ri + R2)/2 = (60™, 0,0). This center-of- 
mass location can be imagined as another "center" of 
another coupled 2-dot system (B) where the two impuri- 
ties are mostly found at the well bottoms of the system 
B. That the relative positions of individual impurities 
vary with respect to their center-of-mass position can be 
imagined as the well barrier height is tuned resulting in 
varying inter-dot separation in system B. A and B only 
interact via the e-I Coulomb interaction. Now, we study 
the influence of the presence of system B on the qubits of 
the system A. We plot in Fig. [51] the singlet-triplet spin 
splitting as a function of one impurity position for sys- 
tem A. We find that J decreases as the impurity which 
stays closer to the system A is positioned further from 
the system A. This effect is expected to be significantly 
enhanced if one applies an external field to engineer the 
impurities in and out the active region of the target sys- 
tem A. As a result, one can see such anti-crossings as 
seen earlier for Z = =f1 where the triplet state stays in- 
between the two lowest singlets. 



0.218 - 




3.5 4.0 4.5 



FIG. 22: (Color online) Singlet-triplet spin splitting J as a 
function of the impurity position in case the two impurities 
(blue circles) are located relatively far from the two potential 
minima, Ri+R2/2=5am and the closest impurity to the ac- 
tive coupled 2-dot system (the grey box) starts at R\ = 3.5am. 
By fixing the center-of-mass coordinate of the two impurities 
at Kc ~ 5am, we imply such two impurities can be considered 
as another coupled 2-dot system which has no tunneling cou- 
pling with the considered 2-dot system in the grey box. That 
the relative positions of each impurity Ri,2 vary with respect 
to the origin (0,0) of the active system can be mapped to 
the problem of other 2-dot system (in open dotted box) with 
varying potential barrier height. Vt of the active system is 
taken to be 30 meV. J is slightly different from that of the 
system without impurity. 
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V. DISCUSSION AND CONCLUSION 

The coupling of two qubits is sensitively affected by 
the presence of a charged impurity. A charged impurity 
which is found inside the right- or the left-dot "breaks 
down" the equality between the right and the left doubly 
occupied singlet states. In case the impurity weakly cou- 
ples with the electrons, the maximally doubly-occupied 
singlet of the dot containing the impurity has a smaller 
energy due to the e-I coupling. This can be considered as 
a tool to distinguish the doped and un-doped QD com- 
ponents in the coupled QD system. 

A relatively strongly repulsive impurity which almost 
equally couples to the two individual QD electrons will 
result in an entire triplet state at B=0 due to the small 
competition between the two e-I couplings. This can be 
observed by positioning the impurity around the origin 
of the double-dot system. When the impurity is found 
close to the bottoms of the confining potential well the 
system will stay in the singlet state where the two elec- 
trons perform the favored singlet state with spins anti- 
parallel to each other. Due to this fact, the influence 
of the inter-dot separation in the presence of impurity 
can be examined by modulating the well barrier VJ, with- 
out changing the single-particle properties in individual 
QDs. The reason is that the relative position of the im- 
purity with respect to the minima of the 2-dot system is 
effectively controlled by changing Vfc. We found that an 
attractive impurity serves to increase the "tunnel" rate 
between the two coupled dots. 

In general, charged impurities "destroy" the maximally 
entangled singlet state by mixing the different singlets 
having varying contributions to the total wave function. 
A strongly perturbative impurity really "messes up" the 
quantum information obtained from this coupled QDs 
because the triplet state "couples" to the second excited 
state by the presence of anti-crossings with a slightly 
larger than zero energy gap. The ground state now can 
favor a doubly occupied state. Because the triplet is al- 
ways obtained with probability 1, it occurs in between 
the two lowest singlets: the exchange singlet tpi and the 
doubly-occupied ^'2 V'f depending on the sign of the 
effective charge Zi. The energy spectrum appears with 
anti-crossings as seen in our results. 

In the presence of two identical charged impurities, the 
breakdown between the two doubly-occupied singlets of 
the left and the right dots is "mended" if the impurities 
are located exactly symmetrically along the line connect- 
ing the two confining potential minima, i.e. the a;-axis. 
As a result, one no longer obtains such a bias between the 
two doubly-occupied singlets in the left- or the right-dot 
as seen in the energy spectrum. 

The mixing of different singlet states and anti-crossing 
coupling of the singlet and triplet states due to the im- 
purity presence results in significant affect on the qubit 
operations, e.g. exchange and C-NOT, in the considered 
coupled 2-quantum-dot system. 

Finally, we discuss the implications of our results for 



the operation of semiconductor spin qubits. Obviously, 
without some quantitative knowledge of impurity loca- 
tions near the qubits, no real comparison between our 
results and experiments would be possible. But we can 
make some general remarks based on statistical consid- 
erations. 

First, semiconductor qubits are typically fabricated by 
lithographic techniques creating QDs from parent 2DEG 
systems. The low-temperature mobility of the 2DEG is 
controlled entirely by the background charged impurity 
density both in Si"^** and in GaAs'^^ systems. Theoretical 
calculations^!^ can provide quantitatively accurate in- 
formation about the ensemble averaged impurity density 
in the 2DEG from the mobility measurements. This in- 
formation about the background impurity density can be 
converted to a statistical probability of finding impurities 
located in specific QD structures since the effective active 
areas of the dots would be known from the lithographic 
structures. For example, typical GaAs 2DEG would have 
10^" to 10^^ charged impurities per sq. cm. whereas Si 
systems are typically dirtier with 10^^ to 10^^ charged 
impurities per sq. cm. This crudely translates to one 
charged impurity every 10 to 100 nm in linear distance 
statistically. Typical lithographic qubits are 20-100 nm 
squares, indicating the presence of 1-5 charged impuri- 
ties per qubit statistically with GaAs (Si) being near the 
lower (higher number) number. Of course, the details of 
impurity locations matter very much, and the statisti- 
cal considerations could be improved and combined with 
our exchange splitting calculations to obtain the effective 
sample "yield" for a given 2DEG mobility, i.e. estimate 
the fraction of samples that would statistically have the 
impurities far enough for spin qubit operations to work. 
But even these elementary considerations show that the 
effective yield is likely to be very low with probably only 
about 1-5% of QD systems being "lucky enough" to have 
the impurities located far enough from the dots for them 
to work as spin qubits. The rest, even before taking into 
account extraneous experimental factors (e.g. leakage, 
noise), simply would not work because there is no effec- 
tive singlet-triplet coupling providing the necessary en- 
tanglement in the system. A possible way of circumvent- 
ing this problem may perhaps be having several electrons 
per dot so that the external charged impurity potential is 
effectively screened by the QD electrons, but this raises 
the problem of having rather weak exchange coupling in 
multielectron dots—. The complex interplay of multielec- 
tron dots in the presence of random charged impurities 
in the background is left as an important future open 
problem in this subject, which may very well be the next 
important step in this direction. 

Second, an important effect of background charged im- 
purities is their fluctuations which produce various charge 
noise signals in semiconductor devices. This charge noise 
may very well be the limiting decoherence mechanism in 
currently existing semiconductor QD spin qubits^. Nor- 
mally, of course, charge noise would not adversely affect 
spin configurations, but semiconductor spin qubits de- 
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pend on the electrostatic exchange coupUng which is af- 
fected by external charged impurities as discussed exten- 
sively in this paper. Therefore, any impurity fluctuations 
would lead to spin qubit decoherence through the charge 
noise mechanism"^^-'^*. Our current work provides quan- 
titative estimate of the strength of such charge noise if 
the fluctuation spectrum (or the fluctuation timescale) is 
known. Deriving microscopic information about charge 
noise from our results could be another future interesting 
direction of research. 

Third, the inevitable presence of static charged impu- 
rities in the background makes every semiconductor dot 
spin qubit, whether in GaAs or in Si (or some other mate- 
rial) , unique since the microscopic electrostatic potential 
environment for each qubit will necessarily differ in a ran- 
dom manner from qubit to qubit due to background im- 
purities as shown in this paper. In particular, the singlet- 
triplet energy separation and the consequent exchange 
coupling will thus be somewhat random in a collection 
of many qubits. This would in turn necessitate charac- 
terization of each qubit in the eventual quantum com- 
puter rather precisely since the gate operations depend 
very strongly on the knowledge of the exchange coupling 
between the dots. Such a characterization will adversely 
affect the scalability of semiconductor spin quantum com- 
puter architectures. Our work shows the importance 
of having ultraclean impurity-free environment even for 
solid state quantum computation, similar to the require- 
ments for topological quantum computation using non- 
Abelian fractional quantum Hall states^, not because 
spin quantum computation needs ultra-high mobilities^, 
but because it requires stable values of exchange cou- 
pling without large qubit-to-qubit variations. Our cur- 
rent work indicates that an order of magnitude reduc- 
tion in the background impurity concentration in GaAs, 
bringing it down to around 10^^ per cubic cm, which 
is also the goal for topological quantum computation^^, 
should lead to the production of 100-1000 spin qubits 
with impurity-free environment rather easily. Unfortu- 
nately, for Si-Si02 systems this is a very stringent con- 
dition because of the invariable presence of large oxide 
charges near the interface^^, but for GaAs^i^ and Si-Ge 
based quantum dot systems^i^i^ low-impurity materi- 
als may become available for quantum computer archi- 
tectures in the near future. 

Fourth, most of the spin qubit manipulation experi- 
ments in semiconductor quantum dots are carried out 
under transport situations using an external dc voltage 
bias. Such external voltage affects the two dots differen- 
tially somewhat similar to the impurity effects discussed 
in the current work since the confinement potentials for 
the two dots are affected differently by the external volt- 
age. Our technique can be used to study this effect. 
Another possible future direction of study could be im- 
purity effects on the full coupled quantum dot energy 
spectra (i.e. in addition to just the double-dot exchange 
energy mainly considered in the current work) since in 
some experiments^ higher energy levels come into play. 



It will, in fact, be interesting also to consider multielec- 
tron quantum dot systems with more than one electron 
per dot and ask how impurity disorder affects the en- 
ergy spectra. This would necessitate a calculation of the 
quantum dot electronic structure using many orbitals per 
dot (e.g. s, p, d, f, ...), which is beyond the scope of the 
current work. 
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Appendix A: Coulomb matrix elements 

1. Single QDs 

To evaluate the e-I Coulomb matrix element, we cal- 
culate the integral: 

where (pn,i{^) is defined in Eq. ([7]). The denominator in 
the right hand side of Eq. (jAip is eliminated using the 
Gaussian identity: 

1^2 r ^-u-r-^^^ (A2) 

r TT Jo 

Integrating integral (|A1[) over 6 and set = t, we arrive 
at: 

-^'.^^./j^ ^;;|;t+^.)! ^^^^ 

Jo 

Next, we use one of the properties of the Laguerre 
polynomials^: 

/ e-'e-P'Ll{at)Li{ht)dt = 
Jo 

r {a) (A + l)„,(/3 + l)„p^" y. (-m),(a), / aV 
mini (A + l)jj! \p) 

{-n)k{3 + a)k fb\'' 
(/3 + l)fefc! \p) 

and the gamma functions for integers to calculate the 
integral in the right hand side of (jA3|) . Finally, we arrive 
at the formula presented in Eq. 

2. Coupled QDs 

In general, the Coulomb matrix elements between the 
electrons and impurity are obtained numerically. Except 
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for few cases with applied constraint conditions to the 
impurity position, it is impossible to obtain a closed al- 
gebraic form for this type of Coulomb interaction. Using 
a similar method as before to calculate the e-I coupling 
elements, i.e. (^I^^-Impl^)' obtain: 



(¥'L(i^)l^e-Impl'^J^(i)) 
for R = (0, 0, z) and 

{'Pl(r)\Vq_i^^\^l(r)) 
(^L(fl)|Ve_implmi)) 



zv^ 



Zyo<^V^e^o-'^«erfc(zo) 



(A5) 



j-^j(xo±ao)^+^j 



.2 «+yt'> ^ 



lo 



for R = (xo, 2/0:0) where oq 
y/lo, and zq = z/Iq. 



am/lo, xo = x/lo,yo 



It is of interest to consider the special situation when 
the impurity is located along the line connecting the two 
potential well minima, i.e. along the x-axis in this case. 
In the above sections, we discussed this case in depth. 
The Coulomb matrix elements for the left and the right 
QD electron are as follows: 



{V>L(R)\VQ.luip\'PLiR)) = ZV^ y/^e 



lo 



(xo ± a^y 



(A6) 



{'PL{R)\Ve_l^p\ipR(L) 



Appendix B: Heitler-London approximation 

The closed analytical form of the exchange energy be- 
tween the two coupled QDs for a single charged impurity 
located arbitrarily in the xy-plane [R = (xq, J/Oj 0)^o] is 
obtained by solving a basic two-level problem: 



a 



+ e 



7(l + ao6J(l + «06,) 



V^Va^ -V2e^" + V2Io{al)-AZe''°-^^Io{ 



2e + oa, ) 



(Bl) 



+ze 2 _ 



2e" 



-2aoxo-xg-yg + 200^0 + xg "f yg 

^ ^Jol t: 



)]}■ 



Here, ag — a/lo, aoa; — lo/^x, CLOy — lo/ly: CLobx — ^o/hx, Appendix C: Singlet-triplet splitting J for 

and ttoby = lo/lby are dimensionless parameters asso- intermediate impurity effective charges 

ciated with a series of lengths such as a,lx,ly,lbx,hy 
These parameters together with depths Vq and Vf, define 
the size and the shape of the confinement double well. 



Heitler-London approximation does not allow one to 
examine the impurity effect on the coupled qubit oper- 
ations manifested in the mixing of the different singlets 
and "coupling" between the singlet-triplet states as al- 
ready observed and discussed in e.g. Fig. ITSTb) The 
impurity-position dependence of Jhl^ as obtained an- 
alytically in the last three terms in Eq. (jBip , now can be 
studied in Fig. [23] as a function of xq . We see that the 
presence of the charged impurity only plays a role as a 
weak perturbative interaction to the total energy even in 
case the effective e-I coupling is found relatively large (the 
case Zi = —1). Such a property like the triplet-singlet 
transition for weak repulsive e-I exchange interaction is 
not found in this case. 



We show the results obtained for several different inter- 
mediate Zi for both the negatively and positively charged 
impurity cases in Figs. [Ml [HI and [231 These plots 
provide a clearer and deeper look into the triplet-singlet 
transition (Figs. [25] and I26p for a strongly repulsive im- 
purity as well as the physics of anti-crossings found in 
the energy spectrum of the coupled 2-dot system. For 
the Zi < case, Fig. [24] serves to explain in detail the 
transition from having a minimum to a maximum in the 
singlet-triplet splitting J between the two lowest energy 
levels. 



23 




0.05 




0.10 



0.05 - 







2 



FIG. 23: (Color online) Singlet-triplet spin splitting J as 
a function of the impurity position along the a;-axis for (a) 
Z\=-Q.l and (b) 0.1. Black solid curves are extracted from 
Fig. [5] and blue dashed curves are obtained from Eq. HB1|I . 
We receive qualitative agreement between our numerical cal- 
culations and analytical results where the singlet-triplet split- 
ting energy calculated using these two methods exhibits a (a) 
minimum for Z\=-Q.l and (b) maximum for Z\ =0.1 around 
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FIG. 24: (Color online) Singlet-triplet spin splitting J as a 
function of the impurity position along the a;-axis for different 
Zi < 0. In this plot one can see the gradual change in J from 
having a minimum around Xi « 2lo to having a maximum 
at different xi depending on Zi as Zi is increases from -0.1 
(black solid) to -0.8 (cyan dash-dotted-dot). 
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FIG. 25: (Color online) Singlet-triplet spin splitting J as a 
function of the impurity position along the x-axis for different 
Zi > 0. The maximum in J is found at different xi as Zi 
changes. 




FIG. 26: (Color online) Singlet-triplet spin splitting J as a 
function of the impurity position along the x-axis for a specific 
case of Zi > 0: 0.6. Inset is the low-level energy spectrum 
of the system. As compared to the pronounced anti-crossing 
seen in Fig. llBf b). the Zi = 0.6 case has the spectrum which 
is close to having pronounced anti-crossings as obtained for 
^1 = 1 in Fig. [Mb V 



Appendix D: Spin splitting energy of Si/SiGe 
coupled quantum dots 



'Si 



Si QDs have a relatively large effective mass ml 
2.8mQ^^g. Plus, the effective Rydberg energy in Si QDs, 
as discussed in Sec. II, is higher than the effective Ryd- 
berg energy of GaAs QDs. These bring into the fact that 
the Si QDs have a smaller electron kinetic energy, which 
means that the electrons tend to be more localized as 
compared to the case of GaAs. Thus the tunnel rate is ex- 
pected to be smaller. Our calculations show that for the 
same model of a coupled 2-dot system in which the shape 
and the size and the barrier height are kept unchanged: 
e.g. Vq = —50 meV and Vf, = 30 meV (correspond- 
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ing to AVb = 9.65 meV), Si/SiGe double-dot system 
has huto « 6.67 meV {Iq w 8.48 nm). The spin singlet- 
triplet splitting J between the two lowest energy levels 
in Si/SiGe double dots is found J « 0.003 meV, much 
smaller than the J « 0.204 meV of the GaAs/AlGaAs 
double dots. In experiments, the electron tunneling rate 
in Si/SiGe coupled double QDs is indeed found much 



lower than the tunneling rate in GaAs/AlGaAs coupled 
QDs (see e.g. Ref."*^). In the presence of a charged im- 
purity, e.g. Zi = —0.1 located at the origin, we obtained 
J increased to J « 0.008 meV for Si/SiGe double QDs 
as compared to J f« 0.26 meV for GaAs/AlGaAs double 
QDs. 
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